Anomalous Viscosity of an Expanding Quark-Gluon Plasma 
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We argue that an expanding quark-gluon plasma has an anomalous viscosity, which arises from 
interactions with dynamically generated color fields. We derive an expression for the anomalous 
viscosity in the turbulent plasma domain and apply it to the hydrodynamic expansion phase, when 
the quark-gluon plasma is near equilibrium. The anomalous viscosity dominates over the coUisional 
viscosity for weak coupling and not too late times. This effect may provide an explanation for the 
apparent "nearly perfect" liquidity of the matter produced in nuclear collisions at RHIC without 
the assumption that it is a strongly coupled state. 



Measurements of the anisotropic collective flow of 
hadrons emitted in noncentral collisions of heavy nu- 
clei at the Relativistic Heavy Ion Collider (RHIC) are 
in remarkably good agreement with the predictions of 
ideal relativistic fluid dynamics 1]. In order to describe 
the data, calculations need to assume that the matter 
formed in the nuclear collision reaches thermal equilib- 
rium within a time Ti < 1 fm/c |3| and then expands 
with a very small shear viscosity <C s, where s is the 
entropy density 01 • The comparison between data and 
calculations indicates that the viscosity of the matter 
cannot be much larger than the postulated lower bound 
"Omin = s/Air Q , which is reached in certain strongly cou- 
pled supersymmetric gauge theories 5] . 

This result is nontrivial because the shear viscosity of a 
weakly coupled, perturbative quark-gluon plasma is not 
small. In fact, the perturbative result for the shear vis- 
cosity, in leading logarithmic approximation, is |^ 

where df ~ O(IOO) is a numerically determined constant 
that weakly depends on the number of quark flavors Uf. 
The result (^), as well as the finding that numerical so- 
lutions of the Boltzmann equation exhibit fluid dynam- 
ical behavior only when the cross section between glu- 
ons is artificially increased by a large factor j3| , have in- 
vited speculations that the matter produced at RHIC 
is a strongly coupled quark-gluon plasma (sQGP). The 
possible microscopic structure of such a state is not well 
understood at present &, 9^ 10] . 

Here we present an alternative mechanism that may 
be responsible for a small viscosity of a weakly coupled, 
but expanding quark-gluon plasma. The new mechanism 
is based on the theory of particle transport in turbu- 
lent plasmas 0, ■ Such plasmas are characterized by 
strongly excited random field modes in certain regimes of 
instability, which coherently scatter the charged particles 
and thus reduce the rate of momentum transport. The 
scattering by turbulent fields in electromagnetic plasmas 
is known to greatly increase the energy loss of charged 
particles 0| and reduce the heat conductivity 0, 



and the viscosity [lil ll^ of the plasma. Following Abe 
and Niu |l7l |. we call the contribution from turbulent 
fields to transport coefficients "anomalous." 

The sufficient condition for the spontaneous formation 
of turbulent, partially coherent fields is the presence of 
instabilities in the gauge field due to the presence of 
charged particles. This condition is met in electromag- 
netic plasmas with an anisotropic momentum distribu- 
tion of the charged particles US'], and it is known to be 
satisfied in quark-gluon plasmas with an anisotropic mo- 
mentum distribution of thermal partons |l9l l2Cl |2l| . 

Most of the work exploring the consequences of the 
instabilities ^22,. .23, . .24] has been focused on the early 
stage of the collision, when the momentum distribution 
is highly anisotropic and far from equilibrium. It was 
pointed out that the fields generated by the instabili- 
ties will drive the parton distribution rapidly toward lo- 
cal isotropy and thus into the hydrodynamical regime 
[isj . Here we are concerned with the later stage of 
the reaction, when the matter is nearly equilibrated and 
evolves by hydrodynamical expansion. Because the par- 
tonic plasma expands rapidly, the momentum distribu- 
tion of the partons remains anisotropic even at late times, 
with the size of the anisotropy being proportional to the 
viscosity. 

As we will show, the turbulent plasma fields induce 
an additional, anomalous contribution to the viscosity, 
which wc denote as 77a- This anomalous viscosity de- 
creases with increasing strength of the turbulent fields. 
Since the amplitude of the turbulent fields grows with 
the magnitude of the momentum anisotropy, a large 
anisotropy will lead to a small value of rj^- Because the 
relaxation rates due to different processes are additive, 
the total viscosity is given by 
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(2) 



This equation implies that rj^ dominates the total shear 
viscosity, if it is smaller than ijc- In that limit, the 
anomalous mechanism exhibits a stable equilibrium in 
which the viscosity regulates itself: The anisotropy grows 
with 77, but an increased anisotropy tends to suppress 77A 
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and thus rj tja- We derive the resulting self-consistency 
condition for tja below. 

The fireballs formed in rclativistic heavy ion collisions 
exhibit collective expansion in both, the longitudinal and 
the transverse directions with respect to the beam axis. 
Here we focus on the longitudinal expansion, but our ar- 
guments apply as well to the transverse expansion com- 
ponent. We assume the longitudinal flow profile during 
the hydrodynamic expansion phase of a relativistic heavy 
ion collision to be approximately boost invariant [2^ and 
of the form Uz{z,t) = z/t, where z and t are measured 
from the collision point in the center-of-mass frame. The 
velocity gradient duz/dz leads to an anisotropy in the 
local momentum distribution 0, 



2T, 



To,- 
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where T denotes the temperature, s the entropy density 
of the matter, and r = — is the time in local 
comoving coordinates. For simplicity, we have assumed 
that the equation of state of the matter is that of free 
massless partons, e = iP = 3sT/4. As (jSJ shows, the 
anisotropy is linearly dependent on the viscosity of the 
matter. 

As already mentioned, any anisotropy of the local mo- 
mentum distribution of quasi-thermal partons engenders 
instabilities of soft gluon modes in the momentum regime 
k < gT, whose growth rate increases in proportion to the 
anisotropy |2ll |. In electromagnetic plasmas the growth 
of the instability eventually saturates when the increas- 
ing field modifies the particle distribution in such a way 
that the instability is eliminated. In nonabelian plas- 
mas the nonlinear self-interactions of the gauge field re- 
strict the growth of the unstable modes: this mechanism 
dominates at weak coupling [2^ |2^ [^Ol ■ In the quasi- 
stationary state that is reached in the nonlinear domain, 
energy absorbed from the thermal partons cascades from 
the most unstable gauge field modes into modes of in- 
creasingly shorter wavelengths. The power spectrum of 
this energy cascade has the form P{k) ^ fc"^, analogous 
to the Kolmogorov cascade in a turbulent fluid • 

As already stated, we are here concerned with the ef- 
fect of the unstable field modes on the transport prop- 
erties of the medium when it has reached the collec- 
tive expansion phase. Following the standard Chapman- 
Enskog theory of transport coefficients, we assume that 
the plasma is driven only slightly out of equilibrium by 
the collective flow and that the local phase-space distri- 
bution can be written as 



(4) 



where /o(p) is the local equilibrium distribution, +(— ) 
applies to bosons (fcrmions), and /i has the form 



fiiP,x) 
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|(5ijV • u) is the traceless 



Here Ay (u) = {SliUj 
part of the flow gradient related to the shear viscosity, 
and A parametrizes the strength of the anisotropy. For 
massless particles, A is related to the macroscopic trans- 
port coefficient of shear viscosity by A = 5ry/s. The boost 
invariant longitudinal expansion corresponds to the per- 
turbation 



h{p) 
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In order to explore the response of the plasma to 
this perturbation, we need to determine the influence of 
the saturated field modes on the propagation of thermal 
plasma particles. The relevant transport theory was de- 
veloped by Dupree ,11.,.12J for an electromagnetic plasma 
in the limit of strong turbulence and weak coupling. We 
now generalize this formalism to a nonabelian plasma. 
Our starting point is the Vlasov-Boltzmann equation for 
the phase space distribution of color charges in a 
color-magnetic field B'^: 



d_ 

dt 



.f{f,P,t)^C[f], 



(7) 



where v = p/E is the velocity of a thermal parton with 
momentum p and energy E, F = gQ°-{E"-+v x B"-) is the 
color Lorentz force, and C[f] denotes the collision term. 
We focus here on the effects of the Vlasov term and re- 
fer to the results of Arnold, Moore, and Yaffe JjJ for the 
viscosity due to incoherent collisions. Because transverse 
field modes have the highest growth rates in the linear 
regime and transverse color-electric fields are less effec- 
tive in restoring the particle distribution to isotropy, we 
here concentrate on the effects due to coherent color- 
magnetic field modes. The effects of longitudinal electric 
field modes, which arc strongly excited in the nonlinear 
domain and can also lead to isotropy, will be discussed 
in a forthcoming longer publication |32l |. 

In order to isolate the dissipative effects of the color 
field, one averages the particle trajectories over an en- 
semble of color-magnetic fields. Assuming that (-B") = 
and factorizing higher than second moments of the field 
distribution, one can then show that the ensemble aver- 
aged phase space distribution / satisfies an equation of 
the Fokker-Planck type [3| '■ 



d_ 

dt 



f^C[f] 



with the diffusion tensor 



dt' {Fdf{t'),t')F,m^t)) ■ 



(8) 



(9) 



Dupree's treatment is based on the argument that 
the Fourier components of the coherent field, B{k, LOk) are 
slowly varying functions and the autocorrelation function 
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of the Lorentz force in ^ is determined by the action 
of the magnetic field on the particle trajectories. This 
leads to a self-consistency condition for the mean de- 
viation (Af^) of the particle trajectories from straight 
lines. Because the growth and coherence of the non- 
abelian gauge field (at weak coupling) is not controlled by 
the back reaction of the distorted particle trajectories on 
the gauge field, but by the inherent nonlinearities of the 
gauge field itself , we here adopt a different approach. 
The velocity of propagation of the collective modes of the 
nonabelian plasma is less than the speed of light, while 
the thermal particles move (nearly) at the speed of light. 
The autocorrelation function of the color-magnetic field 
along the path of a particle will thus be controlled by 
the spatial correlation length of the fields created by the 
growth of the unstable modes. Assuming that the cor- 
relation length for the color-magnetic fields is short in 
comparison with the curvature radius of the trajectory 
of a plasma particle, we can then take v(t) = v(t') = v 
out of the average in ^ and are left with the autocorrela- 
tion function of the magnetic field along a typical particle 
trajectory: 



dt' {Bt{t')B%t))^{BtB))r^ 



(10) 



In our case, the color-magnetic fields generated by the 
plasma instability point in a transverse direction with 
respect to the beam. Assuming that the ensemble av- 
erage is diagonal in color and employing the notation 
—ip X Vp, we can write the diffusive term as 



VpD(p)Vp 



2(iV2 - l)E^ 



(i?^)r^(i^))^ (11) 



where = 3 is the number of colors and the index _L 
denotes the components transverse to the beam axis. 

The action of the diffusion operator on / is easily 
evaluated by noting that the perturbation fi{p) has 
quadrupole form. In order to derive the anomalous vis- 
cosity due to the diffusion term, we follow Abe and Niu 
|l7l | and take moments of the drift and diffusion terms 
in © with p1 . Using massless quarks and gluons in the 
momentum integrals, we obtain: 



d^p 

{2TTfE 



pIv' 
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pl^pD{p)^J{p) 

1 Ag2(i32)r^4iV,^C(2)^2 
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T\ (13) 



where 



Tt (iV2-l) 

16+ 12(1 - 2"^)n/, (14) 
16 + 6(l-2-^)n/(A^2_^)/^2^ (^5) 



and we have used that Q"^ = for gluons and Q"^ = 
{N^ — l)/{2Nc) for quarks. Equating the two results and 
using the relation A — 51] /s we obtain the sought after 
expression for the anomalous shear viscosity due to the 
coherent color-magnetic fields: 



4(7V2 - l)i/4C(4) sT^ 
5iV,i.^C(2) .g2(i?2)T„ 



(16) 



It is noteworthy that the rhs. of (|16|l itself depends im- 
plicitly on the viscosity, because the intensity of the tur- 
bulent fields grows with increasing anisotropy of the mo- 
mentum distribution in the plasma. 

Next we need to address the question how large (B^) 
and are. The coherent color magnetic fields are only 
generated by the plasma instability when the momentum 
distribution of partons in the quark-gluon plasma is de- 
formed due to the collective expansion. We know from 
analytical studies how the growth rate of the instability 
depends on the anisotropy of the momentum distribu- 
tion, but there are no published systematic studies that 
show how the saturation level of the coherent field energy 
depends on the anisotropy. 

The study by Romatschke and Strickland jl^l ex- 
presses the anisotropy in terms of a parameter ^ and 
a unit vector n. Choosing n — Sz, this ansatz corre- 
sponds to a perturbation of the equilibrium distribution 
of the form © with A ~ ^Tt/2. The average intensity 
of the coherent color-magnetic fields is a function of the 
momentum anisotropy. For lack of a precise knowledge 
of this function, we here parametrize it as a power law: 
g^(B^) = bog'^T'^S^^^, and conjecture a linear dependence 
(6i — 1). We argued above that the scale for the mem- 
ory time Tin for the nonabelian plasma will be set by 
the spatial coherence length of the coherent fields. This 
coherence length is given by the wavelength of the maxi- 
mally unstable mode, which is of the order of the Debye 
length: Tm ~ /i^^ ^ (gT)^^ , where the last form assumes 
weak coupling. We currently lack a precise determination 
of Tm, which could be obtained from numerical solutions 
of the Yang-Mills equations for the anisotropy discussed 
here. 

Using the relations between ^, A, and ry, we can now 
state the dependence of the rhs. of p6|l on the viscosity: 



lObog'^. 

ST 



(17) 



Note that t] here is the viscosity due to all sources of 
dissipation, including collisions. If we neglect the viscos- 
ity rjc due to particle collisions and set = t/a on the 
right-hand side of H17|) . eq. (|16|l yields a self-consistency 
condition for the anomalous viscosity ?7a, which has the 
solution: 



2{Nl - l)^4C(4)rT 
256o7VcZ/^C(2) 



1/2 



4 



Several things are notable about this result. First, if 
the memory time Tm is longer than our estimate, the 
value of 77a decreases. Second, the dependence of 77A on 
the gauge coupling is parametrically much weaker than 
that of the coUisional viscosity Thus, for weak cou- 
pling (7 <C 1 and not too late times t, the anomalous 
viscosity will be much smaller than the coUisional viscos- 
ity. According to this implies that rjA is dominant 
at early times, and the coUisional viscosity rjc may dom- 
inate at large times. The cross-over time Tc between the 
two regimes is given by the condition rjA = rjc- Because 
of the failure of the perturbative result for r]c in the do- 
main g > 1, it is difficult to give a reliable estimate of the 
cross-over time between the two regimes. Ignoring those 
limitations and assuming bo ^ 0(1), one can surmise 
that the anomalous viscosity dominate for a few fm/c in 
heavy-ion collisions at RHIC or LHC energies. At late 
times, the transverse expansion of the medium also needs 
to be taken into account. Assuming a radial dependence 
of the form Ur{r) « f3ar/R, where R is the nuclear radius, 
the anisotropy due to transverse flow is 
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p,p,A,,(^.) = -||(3p2_p2) 



(19) 



Comparing with ^ one sees that the effects from the 
radial expansion become comparable to those from the 
longitudinal expansion when /3o ~ R/t. The anomalous 
contribution to the viscosity may, therefore, never be neg- 
ligible during the life-time of the plasma phase. We note 
that collisions among thermal particles may suppress the 
Weibel instability when the coupling constant exceeds 
a threshold value [33l |. It is presently unknown whether 
this occurs before or after the cross-over between rj^ and 
77c for experimentally relevant conditions. 

The approach outlined here can be used to derive other 
anomalous transport properties of an expanding, turbu- 
lent quark-gluon plasma. Maybe the most important 
among these are the coefficient q of radiative energy loss 
of an energetic parton [s^, which might be increased by 
scattering on turbulent fields. 

In summary, we have shown that an expanding quark- 
gluon plasma acquires an anomalous viscosity due to 
the interaction of thermal partons with chromomagnetic 
fields generated by instabilities of soft field modes. In 
the weak coupling limit, the anomalous viscosity is much 
smaller than the viscosity due to collisions among ther- 
mal partons. By reducing the shear viscosity of a weakly 
coupled, but expanding quark-gluon plasma, this mecha- 
nism could possibly explain the observations of the RHIC 
experiments without the assumption of a strongly cou- 
pled plasma state. A definitive answer will require the 
numerical evaluation of the correlation function (|10|l as 
a function of the anisotropy parameter A. 
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